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PACS 76 . 30 . Mi - Color centers and other defects 
PACS 71.70. Jp - Nuclear states and interactions 
PACS 76 . 70 . Hb - Optically detected magnetic resonance (ODMR) 

Abstract - In this Letter, we proposed and experimentally demonstrated a method to detect 
vector magnetic field with a single nitrogen vacancy (NV) center in diamond. The magnetic field 
in parallel with the axis of the NV center can be obtained by detecting the electron Zeeman shift, 
while the Larmor precession of an ancillary nuclear spin close to the NV center can be used to 
measure the field perpendicular to the axis. Experimentally, both the Zeeman shift and Larmor 
precession can be measured through the fluorescence from the NV center. By applying additional 
calibrated magnetic fields, complete information of the vector magnetic field can be achieved 
with such a method. This vector magnetic field detection method is insensitive to temperature 
fluctuation and it can be applied to nanoscale magnetic measurement. 
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' Introduction. — The magnetic field plays a signifi- 
cant role in various regions such as material science and 
quantum physics since it can affect the motions of charged 
particles and the dynamics of electrons and nuclei spin 
states. Therefore, to precisely detect the magnetic field 
is becoming more and more important. Various magnetic 
detection techniques have been developed based on differ- 
ent systems including scanning Hall probe microscope [I] , 
ihagnetic resonance force microscope [2], optical magne- 
tometry [3] and scanning superconducting quantum inter- 
ference device (SQUID) [i]. Very recently, nitrogen va- 
cancy (NV) centers in diamond were proposed to detect 
magnetic field with nanoscale spatial resolution [Ml] ■ The 
general idea of NV magnetometer is that the electron spin 
energy levels change with external magnetic field due to 
electron Zeeman effect. The great advantages of such a 
magnetometer are that the electron spin states in NV cen- 
ter possess long coherence time at room temperature and 
they can be optically initialized and detected jH]. There- 
fore, the NV magnetometry is very potential for future 
practical applications. It is not only convenient to be op- 
erated without cryostat as required in SQUID, but also 
very sensitive to weak magnetic field, even the magnetic 
field produced by single remote nuclear spin [5HTT]. 

It is also important and necessary to detect the vec- 
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tor magnetic field with a single NV center. The interest 
in the detection of vector magnetic field has been driven 
by numerous applications. For example, it can provide 
both translation and rotation information of nanoparti- 
cle in biological cell [T21[T3]. However, the application of 
NV center in magnetic measurement is restricted by the 
fact that the electron Zeeman shift mainly depends on the 
external magnetic field component along the defect's sym- 
metry axis, especially for weak magnetic field. Though 
several NV centers with different orientations have been 
used to detect the local magnetic field vector [T^, they 
were separated by hundreds of nanometers to be optically 
distinguished, which will reduce the spatial resolution. 

Here, we proposed a method to detect the vector mag- 
netic field by a single NV center, with the assistance of ^^N 
nuclear spin. Since the distance between nuclear spin and 
electron spin is less than 1 nm, high spatial resolution can 
be achieved. The magnetic field in parallel with the axis of 
^■^NV center is obtained by detecting the electron Zeeman 
shift. To get the magnetic field perpendicular to the axis, 
the spin-selective pulse is used to measure nuclear spin 
state, which undergoes Larmor precession in the magnetic 
field. The results of electron spin Zeeman effect and nu- 
clear spin Larmor precession can be detected through the 
fluorescence from the ^^NV center. Therefore, the value of 
magnetic field is determined, so is the angle between NV 
center symmetry axis and the magnetic field. Complete 
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Fig. 1: (Color online) (a) The schematic of experimental setup 
of single NV center magnetometry in diamond, (b) Level struc- 
ture of ^^NV center, (c) The ODMR spectrum of ^'^NV center 
with an unknown magnetic field. The spectrum reveals the 
hyperfine structure of ^'''N. 



vector magnetic field can be determined by introducing 
calibrated magnetic fields. Such a vector magnetic field 
sensor can be used for future studies of nano-matcrial and 
biological system. 

Model. — The NV center in diamond consists of a sub- 
stitutional nitrogen atom and an adjacent vacancy with 
Csv symmetry. The ground state of NV electron spin is a 
triplet with 5=1. The Hamiltonian of NV center electron 
spin reads 

H,^ DS^,+gePeB-S. (1) 

Here, z denotes the NV center symmetry axis. S is the 
electron spin vector operator, ge and /3e are the electron g 
factor and Bohr magneton, respectively. D w 2.87GHz is 
zero field split (ZFS) between sublevels nis = and = 
±1 at room temperature (FiglUb)). These spin states 
can be manipulated by microwave and detected through 
fluorescence intensity. Due to Zeeman effect, an external 
static magnetic field will cause the shifts of electron spin 
states. Usually, without the knowledge of the direction of 
external magnetic field, we can not distinguish the electron 
spin transitions between Anis = +1 and Am^ = — 1 only 
from the resonant spectrum. But we can definitely get two 
different resonant frequencies ijj± with a;+ > cj— . Under 
weak field condition with |B| <C D/gePe, the shifts caused 
by magnetic can be expressed as 



Auj± « ±gel3e\B, 



+ n9ePeB±- 



D^-glPl\B, 



(2) 



where B±^ = y -B^ + B'^. Obviously, the second term is 
much smaller than the first term. So the Zeeman shift is 



mainly determined by Bz- 

However, the electron spin of NV center is not isolated in 
diamond. Usually three species of nuclear spins are found 
to be coupled with the electron spins of NV center: ^"'N or 
^^N atom associated with the NV center, and randomly 
placed ^^C atoms. When the distance between nuclear 
spin and electron spin of NV center is small, the nuclear 
spin will affect the dynamics of electron spin significantly. 
The Hamiltonian of a nearby nuclear spin and its hyperfine 
interaction with the electron spin arc 



Hn — — ffriAiB • I 

iJ, = S • A • I. 



I P I, 



(3) 
(4) 



Here, g„ and /3„ are the nuclear g factor and magneton 
respectively. A is the hyperfine interaction tensor. And 
P is the nuclear quadrupole interaction tensor for / > 1 . 

Our magnetic measurement idea is that B can affect 
the nuclear spin indirectly through the hyperfine interac- 
tion, then the dynamics of nuclear spin is more magnetic 
sensitive than the nuclear Zeeman effect. Here, a single 
i5]\fY -^as chosen to present the measurement method. 
The ^^NV center can be found in natural diamond. It also 
can be created artificially by ion implantation |15| . which 
makes it feasible for future application. In such a compos- 
ite system, the Larmor precession of ^^N is affected not 
only by the external magnetic field, but also the electron 
spin state [T51[T7] . Here, we focus on the nuclear free pre- 
cession with Tils = 0. The Larmor frequency is expressed 
as 

WL = |g„,0„Be//| (5) 

with Be/y ~ B • go. The go tensor can be calculated by 
second order perturbation theory [TBllTSl 



go 




A 



yz 





(6) 

where Aij is the element of hyperfine interaction tensor A. 
As the nuclear spin ^^N lies on NV center's symmetry axis, 
the interaction tensor has diagonal form [T7] Aij = AijSij. 
Therefore the relation between nuclear spin Larmor fre- 
quency and external magnetic field can be simplified as 



l.9nAi 



Bh 



(7) 



with the coefficient 



1 



2g./3e 



X Ai 



These measurable frequencies Aw± and w l are functions 
of \Bz\ and \B±_\. Therefore we can solve the absolute 
values of total magnetic field and polar angle between the 
vector and the NV axis. The nuclear Zeeman effect is 
ignored in Aw± since it is much smaller than the electron 
Zeeman effect as Pel Pn ~ 1836. 

Experiment and results. — The experiment was 
performed on a (100) CVD diamond plate from Element 
Six, with nitrogen concentration lower than 5 ppb. The ex- 
perimental setup is shown in Fig|IJa). All measurements 
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Fig. 2: (Color online) (a) Pulse sequence in the experiment, 
(b) ^^N nuclear spin precession. Pink dots are the experimental 
results, and the blue solid line is cosine function fitting. 



were carried out at room temperature with a home built 
optical confocal scanning system. Microwave was trans- 
mitted through a 15 /im diameter copper wire, and a per- 
manent magnet was placed near the sample to provide 
external static magnetic field. 

As shown in Figjljc) , the transition frequencies between 
spin states can be measured through optically detected 
magnetic resonance (ODMR) [51[T1]. When the applied 
microwave was resonant with one of the transition fre- 
quencies, the fluorescence intensity would show a decrease. 
More precise results can be obtained by other methods 
like Ramsey fringes [50]. The ^^N atom with nuclear spin 
1=1/2 split the electron spin states, as shown in FigHJb). 
Therefore, we can observe two valleys in Figlljc) for one 
electron spin transition branch (ms = — > m,, = — 1 or 
rris = — > uis = The Lorcntzian fitting revealed 

two resonant frequencies 2.78465(2) GHz and 2.78774(2) 
GHz, corresponding to different nuclear spin states. Com- 
paring with the zero magnetic field resonant frequencies, 
we got Aw_ = —85.26{6)MHz under the external static 
magnetic field. 

The free precession signal of ^^N nuclear spin was ob- 
tained by two nuclear spin selective microwave tt pulses. 
The pulse sequences and results are shown in Figl^Ja) and 
(b). The electron spin of NV center was first polarized 
to TOs = by 3/is green laser pulses. Without nuclear 
spin polarization, the nuclear spin was assumed to be 
at thermal equilibrium state with the state density ma- 
trix P^\0,+^){0,+^\+Pi\0,-^){0,-^\ following the ba- 
sis \ms,mi){ms,mi\. After the first selective tt pulse (e.g. 
on the transition |0, +^) — > | — 1, +-^)), the density matrix 
became P^\ - l,+i)(-l,+i| + P;|0, -i)(0, During 
the free precession time t, the nuclear spin with rris = 
periodically changed the direction, while the split between 
nuclear spin states in the rus = — 1 subspace prevented the 
precession of nuclear spin [T71HT] . The second microwave 
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Fig. 3: (Color online) The magnetic field intensity imaging 
detected by a single ^^NV. (a) and (b) are the absolute values of 
_B, and Bj^ , respectively, obtained through the electron Zeeman 
effect and ^^N Larmor precession. 



TT pulse mapped the information of nuclear spin onto the 
electron spin. At last, another green laser pulse was used 
to read out the electron spin. The detected photolumi- 
nescence intensity oscillated as I{t) (x Picosi/liruLt) with 
Larmor frequency ujl- Fitting the precession signal by 
cosine function, we got ujl =163.2(3) KHz in Fig. HJb). 

For ^^N. the nuclear g factor is about —0.566. The hy- 
pcrfinc interaction parameters are Axx = Ayy = 3.65(3) 
MHz and A^^ = 3.03(3) MHz [12, thus = Uy = 
— 15.5(1). Substituting these parameters into Eqs. ^ 
and ([7]) , we solved the absolute magnitude of the magnetic 
field as \B,\ = 3.129mr and = 2.426mT. The second 
term in Eq. (2) was about 2AAIHz. It was small but can 
not be neglected for high accuracy magnetometry. Uti- 
lizing this method, we presented the spatial distribution 
of vector magnetic field by moving the permanent magnet 
along X or Y axis (the laboratory coordinate frame) . The 
results of and B± arc shown in Fig. |21 As the angle 
between magnetic field and NV axis varies with position, 
the change of Bz might be contrary to B± at the some 
locations. 

For comparison, the electron spin echo was also mea- 
sured to estimate |B|. As shown in Fig. 4, the spin echo 
signal is modulated by both "'^^N spin and ^"^C spin bath. 
The coupling between ^^C nuclear spin ensembles and NV 
electron spin contributes to the collapse and revival of the 
electron spin echo [TBll^ . while the coupling between ^^N 
and -'^^NV electron spin contributes to the rapid oscilla- 
tions. The result was fitted by the model including the 
hyperfine interaction with ^^N as [H] S* = i-hi(e-(^/^'=)''-f 
6xe-«^-^-)/^=)')x[l-axsm(27rwiT/2)2sm(27rw2T-/2)2], 
where uii and uj2 were Larmor frequencies of "'^^N nuclear 
spin with different electron spin states. The first revival 
peak was at Tre = 24.27(6)/is, corresponding to the revival 
frequency ujrevivai ~ 41.2A"iJz, which equaled to Larmor 
frequency of ^^C. It corresponded to a magnitude of the 
field |B| = iOrevivai I g^^cP^'^c ~ 3.85mT. This agreed with 
the previous results, where |B| 



Bj^ w 3.96mT. 



Combining with electron Zeeman effect, the electron spin 
echo can also provide the information of vector magnetic 
field. However, the detection of spin echo is usually slower 
than that of nuclear free precession. Furthermore, ^'^C nu- 
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Fig. 4: (Color online) Experimental results of electron spin 
echo signal (pink points) and fitting with theoretical formula 
(blue line). Inset: the experimental control pulse for electron 
spin echo signal. 

clear spins are randomly distributed around the NV cen- 
ter, and the nuclear spins close to the NV center will lead 
to very rapid oscillations in the spin echo signal. All of 
these make the revival peaks difficult to be located pre- 
cisely and lead to detection error. 

Discussion. — Now we will discuss the sensitivity in 
Figs[T]and[5] The ODMR spectrum has Lorentzian shape, 
and there are two points on the spectrum with the highest 
sensitivity. The noise of resonant frequencies at these two 
points are GOKHz. Therefore the noise of Bz caused by 
ODMR method is about 2fiT. 

To analyze the sensitivity of B±, Eq.Q is rewritten as 

ax v \gnPnr 

In this experiment the noise caused by nuclear spin pre- 
cession is much larger, so only the error of a;^ is taken into 
account here. And the precession signal is well fitted by 

I{t) = /o + /cCos(27rwLt)e"*/^° , (9) 

where Tq = 156/is is the nuclear spin coherence time in the 
present measurement. The maximal sensitivity is obtained 
at cos(2niL)Lt) = The uncertainty SB± is 

SB^ = SI^ « 61- ^^^e*/^^ (10) 

After polarizing the nuclear spin by selective tt microwave 
pulse and optical excitation pulse (results not shown) [T7] , 
the noise ratio is reduced to 4^ « 0.045. Then the error of 
B± with our method is about 1.8/xT at t = Tq. Therefore, 
the minimal detectable amplitude of magnetic field is at 
scale of fiT with our method, and the error of angle is lower 
than 0.05°. The NV center used in our experiment is very 
deep inside the crystal (about decades of micrometers), 
then only 3 x lO'^ photons are detected per second. Using 
the implanted ^^NV centers which are near the sample 
surface, the collection efficiency can be improved. And 
longer coherence time of ^^N nuclear spin will also improve 
the sensitivity. 



With the method in Ref . [S] , where all electron spin res- 
onant frequencies' shifts Auj± are measured with ODMR, 
\Bz\ and \B±\ can also be calculated by Auj+ — Aw_ and 
Aw+-|-Aw_. In this way, the error 6B± « D gg2^2^~ i will 
be about 30/zT with B± = 2AmT. It is higher than our 
method. Further more, the ODMR method is affected by 
the fiuctuation of temperature. If the temperature changes 
1 K at room temperature, D will shift about 84 KHz . 
Then w+ + lj- changes 168 KHz, and the corresponding 
error of _Bx could be more than lOO/iT in weak magnetic 
situation. In contrast, the Larmor frequency of ^^N is 
mainly determined by B±^ , and it is insensitive to temper- 
ature fluctuation. 

Although only the magnitude and polar angle of the 
vector magnetic field are obtained so far, it is possible to 
deduce all the vector information by additional calibrated 
fields. As shown in FiglS] (a), all possible magnetic vec- 
tors started at a fixed point in the vector space with the 
endpoints limited on two rings (in blue). By adding the 
additional calibrated field, we can get another two rings 
of the possible endpoints of the detected magnetic. The 
two yellow rings in FigIS] (a) show the measured results 
by subtracting the additional calibrated field. If the cal- 
ibrated field is neither perpendicular nor parallel to the 
NV axis, there would be only two possible directions of 
the unknown magnetic field left. In this way, two differ- 
ent calibrated fields are needed to calculate the complete 
magnetic field vector. 

We also propose a near field scanning vector magne- 
tometry with nanoscale spatial resolution, as shown in 
Fig. 5(b). A single ^^NV center is fixed at the tip of 
atomic force microscope (AFM) probe. There are metal 
strip lines on the tip which can generate calibrated mag- 
netic field with fixed direction by applying specific cur- 
rent. As the distance between nuclear spin and electron 
spin is lower than 1 nm, the spatial resolution should be 
several nanometers, which is limited by the AFM. Using 
this scanning magnetometry, the three-dimensional vector 
magnetic field distribution can be measured with high pre- 
cision and spatial resolution, which can be used to study 
the dynamics of biological systems, and probe the proper- 
ties of nanostructure which consists of magnetic materials. 

Conclusion. — In summary, we proposed and demon- 
strated a vector magnetic field magnetometry by a single 
NV center with an ancillary nuclear spin. In this compos- 
ite spin system, the frequency shift of electron spin and the 
Larmor frequency of nuclear spin can be detected optically 
when the NV center is placed in external static magnetic 
field. Then, the amplitudes of magnetic field parallel and 
perpendicular to the axis of NV are determined. Further 
more, by applying calibrated magnetic fields with defined 
directions, complete information of the vector can be mea- 
sured. In future work, the precision of measurement can be 
improved by methods like dynamic decoupling. With high 
spatial resolution system like AFM, the magnetic proper- 
ties of nanoscale particle can be detected. 
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Fig. 5: (Color online) (a) The results without calibrated field 
defines two rings (blue) in the vector space, as the azimuthal 
angle is not sure. Another two rings (yellow) can be obtained 
by applying a calibrated field. The intersections of the results 
with and without calibrated field indicate that the real direc- 
tion (orange) is in the two possible directions, (b) A proposal 
schematic of a compact magnetic probe tip. 



* * * 

This work was supported by the National Funda- 
mental Research Program of China under Grant No. 
2011CB921200, the Knowledge Innovation Project of Chi- 
nese Academy of Sciences (60921091), National Natural 
Science Foundation of China under Grant No. 11004184. 



REFERENCES 

[1] Chang A. M., Hallen H. D., Harriott L., Hess 
H. F., Kao H. L., Miller R. E., Wolfe R., van der 
ZiEL J. and Chang T. Y., Appl. Phys. Lett, 61 (1992) 
1974. 

[2] Sidles J. A., Garbini J. L., Bruland K. J., Rugar 
D., ZijGER O., HOEN S. and Yannoni C. S., Rev. Mod. 
Phys., 67 (1995) 249. 

[3] BUDKER D. and Romalis M., Nature Phys., 3 (2007) 227. 

[4] Bending S. J., Adv. Phys., 48 (1999) 449. 

[5] Balasubramanian G., Chan I. Y., Kolesov R., Al- 
Hmoud M., Tisler J., Shin C, Kim C, Wojcik A., 
R.Hemmer p., Krueger A., Hanke T., Leitenstor- 
FER A., Bratschitsch R., Jelezko F. and Wrachtrup 
J., Nature, 455 (2008) 648. 

[6] Taylor J. M., Cappellaro P., Childress L., 
Jiang L., Budker D., Hemmer P. R., Yacoby A., 
Walsworth R. and Lukin M. D., Nature Phys., 4 (2008) 
810. 

[7] Maze J. R., Stanwix P. L., Hodges J. S., Hong S., 
Taylor J. M., P.Capperllaro, L.Jiang, Dutt M., 
ToGAN E., ZiBROv A. S., Yacoby A., Walsworth 
R. L. and Lukin M. D., Nature, 455 (2008) 644. 

[8] Manson N. B., Harrison J. P. and Sellars M. J., 
Phys. Rev. B, 74 (2006) 104303. 

[9] Zhao N., Honert J., Schmid B., Klas M., Isoya J., 
Markham M., Twitchen D., Jelezko F., Liu R.-B., 



Fedder H. and Wrachtrup J., Nature Nanotechnology, 
7 (2012) 657. 

[10] KoLKOwiTZ S., Unterreithmeier Q. p., Bennett 
S. D. and Lukin M. D., Phys. Rev. Lett., 109 (2012) 
137601. 

[11] Taminiau T. H., Wagenaar J. J. T., VAN der Sar T., 
Jelezko F., Dobrovitski V. V. and Hanson R., Phys. 
Rev. Lett, 109 (2012) 137602. 

[12] Hall L. T., Heart G. C. G., Thomas E. A., Simpson 
D. A., McGuiNNESS L. P., Cole J. H., J. H. Manton 
R. E. S., Jelezko F., Wrachtrup J., Petrou S. and 
HOLLENBERG L. C. L., Sci. Rep., 2 (2012) 401. 

[13] McGuiNNESS L. P., Yan Y., Stagey A., Simpson 
D. A., Hall L. T., Maclaurin D., Prawer S., Mul- 
VANEY p., Wrachtrup J., Caruso F., Scholten R. E. 
and HOLLENBERG L. C. L., Nat. Nano., 6 (2011) 358. 

[14] MAERTZ B. J., WlJNHEIJMER A. P., FUCHS G. D., 

NowAKOWSKi M. E. and Awchalom D. D., Appl. Phys. 
Lett, 96 (2010) 092504. 
[15] Rabeau J., Reichart P., Tamanyan G., Jamieson D., 
Prawer S., Jelezko F., Gaebel T., Popa L, Domhan 
M. and Wrachtrup J., Appl. Phys. Lett, 88 (2006) 
023113. 

[16] Childress L., Gurudev Dutt M. V., Taylor J. M., 

Zibrov a. S., Jelezko F., Wrachtrup J., Hemmer 

P. R. and Lukin M. D., Science, 314 (2006) 281. 
[17] Dutt M. V. G., Childress L., Jiang L., Togan E., 

Maze J., Jelezko F., Zibrov A. S., Hemmer P. R. 

and Lukin M. D., Science, 316 (2007) 1312. 
[18] Maurer p. C, Maze J. R., Stanwix P. L., Jiang L., 

GoRSHKOv A. v., Zibrov A. A., Harke B., Hodges 

J. S., Zibrov A. S., Yacoby A., Twitchen D., Hell 

S. W., Walsworth R. L. and Lukin M. D., Nature 

Phys., 6 (2010) 912. 
[19] Cruder A., DrSbenstedt A., Tietz C, Fleury L., 

Wrachtrup J. and C.von B., Science, 276 (1997) 2012. 
[20] NuSRAN N. M., Momeen M. U. and Dutt M. V. G., 

Nature Nanotech., 7 (2012) 109. 
[21] Dreau a.. Maze J.-R., Lesik M., Roch J.-F. and 

Jacques V., Phys. Rev. B, 85 (2012) 134107. 
[22] Felton S., Edmonds A. M., Newton M. E., Mar- 

TiNEAu p. M., Fisher D., Twitchen D. J. and Baker 

J. M., Phys. Rev. B, 79 (2009) 075203. 
[23] Rowan L. G., Hahn E. L. and Mims W. B., Phys. Rev., 

137 (1965) A61. 
[24] Chen X.-D., Dong C.-H., Sun F.-W., Zou C.-L., Cui 

J.-M., Han Z.-F. and Guo G.-C, Appl. Phys. Lett, 99 

(2011) 161903. 



p-5 



